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Summary. We propose a new method for determination of element abundances 
in stellar atmospheres aimed for the automatic processing of high-quality stellar 
spectra. The pan-spectral method is based on weighted cumulative line-widths 
Qx — J* | 1 (1 — R\)d\, where R\ is residual flux and Z is abundance of stud- 
ied element. Difference in quantities Q\ found from synthetic and observed spectra 
gives a correction to the initial abundance. Final abundances are then found by 
rapidly converging iterations. Calculations can be made for many elements simul- 
taneously and do not demand supercomputers. 

1 Description of the method 

Essential developments in observational high-precision and high-resolution 
spectroscopy, and fast-growing computing facilities stimulate to seek new au- 
tomatic methods of analysis of stellar spectra. We recently started to develop 
a pan-spectral (or broad-band) method for determination of abundances of 
chemical elements and their isotopes in stellar atmospheres and to compose 
a software necessary for its application. 

The method is based on weighted cumulative line-widths Q\ defined as 



where R\ is the residual flux (intensity) and Z — log(N e i em /N to t) is the 
abundance of studied element or isotope. 

Integral (1) can be easily reduced to the equivalent width of a spectral line 
by omitting the derivative \dR\/dZ\ and integrating only over wavelengths of 
this line. We call the quantity Q\ cumulative because integral is taken over 
all lines starting from some initial wavelength to the current wavelength. And 
we call it weighted because the derivative of residual flux R\ with respect to 
abundance Z can be treated as a weight function. This derivative automat- 
ically excludes spectral regions insensitive to changes of the abundance of 
studied element and gives a large contribution in the most sensitive regions, 
i. e. in the centers of non-saturated lines and in the steep wings of strong 
lines of this element. 




(1) 
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Weighted cumulative widths Q\ for both observed and synthetic spectra 
are calculated using the derivative of R\ found from the synthetic spectra. 
For small abundance changes the derivative can be replaced by differences, 
thus 

dR\ _ R x - R\ 



dZ Z+ - z- 



(2) 



where flux R^ corresponds to abundance Z + = Z + SZ and flux corre- 
sponds to abundance Z~ = Z — SZ. Correction to the abundance of studied 
element Z can then be found by formula 

where Q°^ s is found from the observed spectrum and Q s ^ n from the syn- 
thetic spectrum with initial (guessed) abundance Z. Quantities Q\ and 
correspond to synthetic spectra with abundances Z + and Z~ , respectively. 

Abundance of element can be found iteratively using the corrections AZ, 
the number of iterations needed depends on difference between initial and 
best-fit abundance values. Abundance determined by such a procedure is a 
best fit for all lines of the element, diminishing thus an influence of possible 
atomic data errors. We would like to point out that here the contribution 
of blended lines is duly taken into account. Abundances can be found simul- 
taneously for many elements. Our calculations show that the method works 
well for line-rich elements, but for line-poor elements further refinements are 
necessary. 

The method can be used also for determination of the effective tempera- 
ture and gravity, using as the weight function the derivative of residual flux 
with respect to T e s or logg, respectively. 



2 Testing the method 
2.1 Observational Data 

We applied the weighted cumulative width method to determine element 
abundances in the atmosphere of the chemically peculiar star HD f75640. 
High resolution spectral atlas of this star, covering the 3 040 - 10 000 A re- 
gion, was recently published by Castclli and Hubrig [1]. The spectrum with 
resolution 90 000 - 110 000 and signal-to- noise ratio 200 - 400 was obtained at 
ESO 8 m UT2 telescope with the UV- Visual Echelle Spectrograph (UVES). 
Castelli and Hubrig [1] also computed a synthetic spectrum for the whole ob- 
served interval with the SYNTHE code [4] and determined abundances for 48 
ions of 39 elements by the visual best fit of synthetic and observational spectra 
of selected lines. The model atmosphere accepted was an ATLAS12 model [2] 
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with parameters T e fj = 12 000 K, logg = 3.95, km s -1 . Improved and ex- 
tended Kurucz atomic line list used in computations is available at F. Castelli 
website (http : / / wwwuser . oat . ts . astro . it/castelli/ stars . html). 

We had to ignore most of the IR spectrum due to telluric lines and also 
some spectral regions with weak overlap or gaps between adjoined echelle 
spectral orders. 

2.2 Synthetic spectrum 

In computations of synthetic spectra we based on the Kurucz ATLAS9 [3] 
model atmospheres and on the line lists and initial element abundances pro- 
vided by Castelli and Hubrig [1]. Computations of the synthetic spectra 
were carried out with program SMART composed by our team during last 
decade [5] . SMART is a compact and simple software and does not demand 
large computer facilities. We used the same spectral resolution (500 000) as 
Castelli and Hubrig in computations with SYNTHE [4]. Gauss broadening 
and the noise cutoff parameters were found to achieve best fit of observed 
and synthetic spectra. Obtained synthetic spectrum agrees reasonably well 
with the observed one. 

2.3 Application of the pan-spectral method 

To estimate the domain of applicability of the method, we obtained abun- 
dances of 12 line-rich elements (Mn, Ti, Cr, Fe, Y, Mg, Si, Ca, O, S, Yb, Ni) 
in the atmosphere of HD 175640. Regions of Balmer and Paschen continua 
were studied separately because abundances of studied metals turned out 
to be larger in outer layers of atmosphere (UV region) than in inner layers 
(visual region). 

We started with abundances derived by Castelli and Hubrig [1] and iter- 
atively found corrections to these values. Abundances of Si, O and S from 
UV spectrum have not been found since there were only few weak lines of 
these elements in UV region. Final abundances were found with only two it- 
erations since necessary corrections were small. Weighted cumulative widths 
Qx for chromium computed with initial abundance are presented in Fig. 1. 
Corresponding corrections to Cr abundance are +0.16 dex in UV and -0.14 
dex in visual region. 

Obtained abundances of the studied elements and comparison with re- 
sults of Castelli and Hubrig [1] arc presented in Table 1. Our results differ 
somewhat from the abundances obtained by Castelli and Hubrig. However, 
the differences in obtained element abundances may partly be caused by use 
of different spectral synthesis programs, not only by difference of used meth- 
ods. This is particularly true for S and Si. Also the difference of used model 
atmospheres should be mentioned. 

We conclude that the proposed method is promising for analysis of chem- 
ical composition and other main parameters of stellar atmospheres. 
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Table 1. Element abundances [N/N®] for HD175640 





Castelli & Hubrig (2004) 


Present 


study 




UV 




Vis 




UV 


Vis 




ion I ion II 


ion III 


ion I ion II 


ion III 






Mn 


+2.47 +2.40 




+2.46 




+2.75 


+2.57 


Ti 


- +1.43 




- +1.30 




+1.62 


+1.31 


Cr 


+1.19 +1.03 




+1.13 +0.96 




+ 1.24 


+0.96 


Fe 


-0.36 




-0.21 -0.30 




-0.23 


-0.33 


Y 


- +3.38 




- +3.01 




+3.26 


+3.04 


Mg 






+0.18 -0.25 




-0.19 


-0.09 


Si 






- -0.23 


-0.09 




+0.21 


Ca 


- -0.15 




+0.42 +0.06 




+0.10 


+0.24 


O 






+0.03 






+0.10 


s 






- -0.41 






+0.03 


Yb 


- +3.14 


+3.66 


- +2.76 




+3.64 


+2.79 


Ni 


- -0.22 




- -0.35 




-0.12 


-0.43 
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